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We present a method for approximating the effective consequence of generic quantum gravity 
corrections to the Wheeler-DeWitt equation. We show that in many cases these corrections can 
produce departures from classical physics at large scales and that this behaviour can be interpreted 
as additional matter components. This opens up the possibility that dark energy (and possible 
dark matter) could be large scale manifestations of quantum gravity corrections to classical general 
relativity. As a specific example we examine the first order corrections to the Wheeler-DeWitt 
equation arising from loop quantum cosmology in the absence of lattice refinement and show how 
the ultimate breakdown in large scale physics occurs. 

PACS numbers: 98.80.-k, 95.36.+X, 95.35. +d, 95.85. Ry 



Quantum cosmology can be studied within the con- 
text of mini-superspace models, reducing the full quan- 
tum field theory to a quantum mechanical system of fi- 
nite degrees of freedom. Applying this to the evolution 
of a three metric results in the Wheeler-DeWitt equation 
(WDW) equation which breaks down near the classi- 
cal big-bang singularity. Loop Quantum Gravity (LQG), 
another approach to canonical quantisation, uses triads 
and connections rather than metrics and extrinsic curva- 
tures as the basic variables, with more success [2J, whilst 
string and brane theories remove the singular behaviour 
by the extra dynamical states that become available at 
small scales. Irrespective of the full underlying theory 
the WDW equation must be recovered as a semi-classical 
approximation to the dynamical equations of the theory, 
if classical general relativity is to be produced at large 
scales. We can then ask whether Quantum Gravity (QG) 
corrections can have any effects on large scale physics. 

At first sight it would appear unlikely that this is pos- 
sible. However, the wave-function solutions to the WDW 
equation oscillate and so their derivatives can become sig- 
nificant at large scales. If the QG corrections depend on 
the derivatives of the wave-function, then it is possible 
that they become significant at macroscopic scales. This 
is the case for Loop Quantum Cosmology (LQC), which 
predicts that the evolution of the universe is dictated by 
a difference, rather than a differential equation. If the 
lattice size of this difference equation is fixed, as was as- 
sumed until recently [H, not all the oscillations of the 
wave-function can be supported, leading to a breakdown 
in large scale classical physics d, HJ ■ 

We investigate the effects of a general class of cor- 
rections to the WDW equation, considering all possi- 
ble derivative terms and show that, at least initially, 
these corrections mimic the behaviour of additional mat- 
ter components. This raises the exciting possibility that 
QG corrections could, at late times, produce a cosmo- 
logical constant like effect. It is also possible that such 
corrections may produce additional matter components 
such as dark matter, however this would only be true for 
a limited class of corrections. What is generally true is 
that QG correction terms that dominate at small scale 
can also produce significant large scale effects, which can 



be well approximated by the classical behaviour of addi- 
tional matter sources. 

By considering the mini-superspace model of an ho- 
mogeneous and isotropic universe, the WDW equation 
reads 
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where a is the scale factor, k = 0, ±1 is the curvature, Tirf, 
is the matter Hamiltonian; in our units H = c = 1 [ll| . 

In general, TL^, will have several terms, with different 
dependence on a and ip. Being only interested in the large 
scale (a 3> ^pi) behaviour of Ti.^, we expect one term 
to dominate and we approximate it by [3/(2Trlp- l )]H c /, — 
e ((f)) a 5 . With this approximation the WDW equation can 
be written in as 
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We only consider the flat case, so that solutions to Eq. ([2]) 
are analytically tractable. If this equation breaks down at 
small scales, due to some QG effects, extra terms should 
become significant on these scales. These new QG correc- 
tion terms will be a function of ip (a) and its derivatives, 
so the full underlying equation can be written as 
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where ao is the scale at which the new terms become 
important. At large scales, a 3> ao, one might naively 
assume that the QG corrections can simply be ignored, 
however this assumes that ip (a) and all its derivatives 
remain small at all scales, which is easily shown not to 
always be consistent. This has been explored within the 
context of LQC @,H@, 

in which the exact form of the 

ao corrections are known. 

Assuming ao <§; e (</>), whilst / (ip (a) , d a ip, ■ ■ ■ ) is 
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TABLE I: The scaling of the energy density p and the Hamil- 
tonian Ti^ — a 3 p for matter, radiation and vacuum energy 
with respect to both a and p. The parameter 8 is also given. 
Note that the energy density and hence of vacuum energy 
are negative. 



small, we can find approximate solutions to Eq. ([3]) as 
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where J and V are Bessel functions of the first kind and 
second kind, respectively, C\ and C 2 are integration con- 
stants; for clarity the <fi dependence has been suppressed. 
We can thus evaluate approximations to the derivatives: 
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etc. 
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Taking the large argument expansion of the Bessel func- 
tions, which is the large a limit for S > —3, one obtains 

Ci sin (A (q))-C 2 cos (A (a)) 
hm Z x (a) = - , (7) 



C\ cos (A (a)) + C 2 sin (A (a)) 
r (a)" 
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where V>* (a) is the solution corresponding to choosing the 
integration factors C\ — > — C 2 and C 2 — > Ci compared to 
the solution given in Eq. Q . This approximation is valid 
only when the correction terms are still small compared 
to the matter component, however this method will give 
the correct effective initial behaviour of QG corrections. 
In particular, the corrections become non-negligible for 
scales at which the corrections are of the order of a^ 1 
and the approximation breaks down when they are of the 
order of e (4>) a s . It is clear that for 5 > —1 the derivative 



terms grow, in agreement with Ref. @- In addition, for 
5 > —1, the terms containing the highest powers of e 
grow fastest and hence are always dominant, implying 
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etc. 



This allows us to approximate the initial effect of large 
scale break down of classical physics (breakdown in pre- 
classicality 0]) as a second matter component. In gen- 
eral the (3 th derivative is proportional to a/ 3 ( 1+l5 " 2 , thus 
correction terms like a a d^^p (a) can be approximated by 
additional terms in Eq. ©, that scale like a a -+l 3 ( 1 + s )/^ 2 . 
It is then possible to consider these correction terms as 
being produced from an effective Hamiltonian. A useful 
pedagogical example is to consider classical matter terms 
(although the method is, of course, general), where the 
effective Hamiltonian that mimics the correction terms is 



Hcor (X a 
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An example of how these approximate solutions compare 
to the exact case is given in Fig. (JTJ. 

We have shown how hypothetical correction terms to 
the WDW equation can mimic the behaviour of addi- 
tional matter sources. In particular, such corrections can 
be used to produce an effective cosmological constant. 
For example, consider Ji^ oc a , i.e., a matter dominated 
universe (see Table[I]); 5 — 0. If we want the quantum cor- 
rections to mimic a vacuum energy then we need 7i C or <x 
a , implying a + (3/2 = 3. Thus, correction terms like 
clqcl 3 , aoa 2 d 2 Vyda 2 , a^ad ip(a)/da , aod 6 ip(a)/da 6 , etc., 
with ao < 0, all resemble a vacuum energy in the pres- 
ence of a matter dominated universe (see Fig. {!])). In 
the presence of other types of matter S changes and so 
the form of the correction terms that give vacuum energy 
like behaviour are different. In particular, for a universe 
dominated by radiation, we have S = — 1 and the only 
correction term that can mimic the behaviour of dark 
energy is aga 3 . Considering a universe dominated by a 
field that scales like H,p oc a 2 , i.e., 5 = 2, then a cor- 
rection term like aQa~ 3 d i ip{a)/da i mimics a (negative) 
cosmological constant like term. This is precisely the 
dominant correction from LQC. 

Up to now we have only discussed how the correction 
terms scale with a, however there is clearly also a depen- 
dence on the constant part of the matter Hamiltonian, pa 
(through e (</>)). In principle this differentiates these cor- 
rection terms from other models of dark energy (cosmo- 
logical constant, quintessence, etc.), as does the fact that 
as these correction terms begin to dominate the WDW 
equation this approximation will break down and the be- 
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FIG. 1: Solutions to the WDW equation, for 5 = (i.e., 
matter dominated universe), e{4>) = 5 x 1(T 3 , a = 9 x 1(T 6 
(with lp\ — 1). The QG corrections are aoa 2 6 2 ip (a) /d 2 a. The 
solutions are calculated numerically: the bare solution (solid) 
excludes QG corrections entirely whilst the approx. solution 
(dots) approximates them as discussed in the text. For small a 
the QG corrections do not play a role, however as a becomes 
larger the effects of the QG corrections become significant. 
The form of this deviation is well characterised by the addition 
of a second matter component, which is demonstrated by the 
accuracy of the approximate solution (see dotted line). For 
a ~ 165, the approximate correction term is « 0.25eip (a) 
and hence can no longer be considered small compared to the 
matter component, eip (a); the approximation breaks down. 



haviour of the wave- function will be drastically altered. 

In addition to corrections that mimic the behaviour of 
dark energy, it is also possible to find correction terms 
that have a dark matter like form, for which Eq. ([10]) is 



H cor oc a a+ ^ 1+ ^ 3 oc a 
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For example, in a matter dominated universe (i.e., 6 = 0), 
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all produce additional matter like terms. Notice that, un- 
like the dark energy case, these corrections do not scale 
faster than original matter component and so will never 
dominate Eq. ([3]). In addition, these correction mat- 
ter terms will be closely related to the physical matter 
Hamiltonian degrees of freedom. For example, 



aoa 



- 2 dV(a) 
da 4 



a o e(0) ip(a) 
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which amount to replacing e (0) — > e (4>) = e (4>)+a e (tfi) 
in Eq. ([3]). Within the mini-superspace model used here 
it is only possible to explore homogeneous, isotropic so- 



lutions and so we cannot say that these dark matter like 
correction terms would produce the necessary behaviour 
to explain structure formation, galaxy and cluster dy- 
namics, etc. Another difficulty with this type of correc- 
tion comes from the fact that the energy density of dark 
matter is approximately five times that of standard mat- 
ter. If correction terms were to produce dark matter in 
the presence of a matter dominated universe, ao would 
have to be fine tuned to ensure that ao <C e (4>) (the ap- 
proximation used here) and aoe (0) 2 1 ~ 5. Similar 
dark matter like correction terms can be produced for a 
radiation, or vacuum energy, dominated universe; if these 
corrections arose due to the the presence of a vacuum en- 
ergy, one may explain the coincidence problem, i.e., that 
the dark matter degrees of freedom would be dictated by 
those of the dark energy, however again significant tuning 
would probably be required. 

LQG is a background independent, non-perturbative 
method of quantising gravity. Reducing the symmetries 
of the theory to a specific cosmological model makes 
the theory tractable and ensures a large scale continuum 
limit. The theory is based on holonomies of the standard 
Ashtekar variables , namely the triad and connections. 
In an isotropic model these can be parameterised by sin- 
gle variables, p and c respectively, which in terms of stan- 
dard cosmological variables are \p\ — a . , c = k + 7a, 
(7 is an ambiguity parameter known as the Barbero- 

Immirzi parameter). Define p = pV^ 3 and c = cVg 1//3 , 
where Vb is the volume of a fiducial cell Q, related 
via the classical identity, {c,p} = 87rZpj7/3. By ana- 
logue with the full LQG theory, c is quantised via its 
holonomies, h = exp(i/zoc/2), where fio is an arbitrary 

real number. Then p\/i) = \p\\(i) = ^^p^lfi), and 

e/W2|^ = e^&lfi) = \fi + no). 

The action of the gravitational part of the Hamiltonian 
constraint on the basis states reads [9|, 



where 



-{[S(/i) + S(/i + 4 Mo )]|m + 4mo} 



2567^^3 
-4%) + [S{ji) + S(fi - 4 M0 )] |M - 4 M0 ) } (14) 



SQ*) = (4^ 2 l7 /3) 3/2 | M + A10 | 3/2 -|A*-Mo| 3/2 • (15) 

The Hamiltonian constraint is (li. g + Ti-^j \4>) = 0, where 

\ip) = VmIm)- Taking the continuum limit — > (m) 
and using the expansion [4| 



| M + a/i | 3/2 -|/x + ^o| 3/2 = /i 3/2 
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we can expand Eq. (fT4]) . 

Changing variables, using a 2 = 47rJpj7|/z|/3, we obtain 
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where ao = lQn 2 l^fi 2 ij^/9. The full Hamiltonian con- 
straint is easily written in the same form as Eq. ©. As 
we have shown, for 6 > — 1 the higher derivative terms 
scale faster. Since in this case the lower derivative terms 
are further suppressed by powers of a, to a good approx- 
imation the effective large scale equation to solve is 




l_d 

a da 



[aip (a)] 
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1 dip (a) 


da 


a da 



+ {<pf a 2S - V (a) + e (0) ofy (a) « , (18) 

where 3H < p/(2nlp l ) = e((p)a s and we used Eq. (Si) to 
approximate the fourth order derivative term in Eq. (117|) 
by term that resembles an additional matter component. 

LQC is a concrete example of the general methods for 
approximating the QG corrections. For a universe dom- 
inated by a (classical) matter content with 5 = 0, the 
leading correction term acts as an effective radiation field, 
T^cor — a o e (0) 2 a_1 /3, i-e., the correction term mim- 
ics radiation. For a radiation dominated universe with 
6 = —1, the correction term acts as TC CO r = Q-o e {<P) 2 a, 3 - 
For 6 = 3, i.e., a universe dominated by a vacuum en- 
ergy, such as during inflation, the correction term acts 
like Hcor — clq£ (4>) ° 5 - Notice that in this case the 
corrections scale faster than the existing matter com- 
ponent, which motivated the modelling of lattice refine- 
ment (fiQ — > ft (u) is no longer a constant) in loop quan- 
tum cosmology 0, S ■ For a universe dominated by a 
source with 6 — 2, the dominant correction term acts as 
Hcor = a o e {4>) fl3 i i- e - h mimics a negative cosmological 
constant (i.e. an accelerated contraction), whilst 6=1/2 
leads to the correction term mimicking dark matter. 

We have shown how generic QG corrections to the 
WDW equation mimic the behaviour of additional mat- 



ter sources, at least whilst they are smaller than the ex- 
isting matter components. This simple procedure can be 
used to examine, to first order, the effect of any QG cor- 
rection under consideration. Here, as a specific example, 
we used the first order corrections from loop quantum 
cosmology to examine what the breakdown of large scale 
classical physics looks like in the absence of lattice re- 
finement. The fact that quantum corrections introduce 
extra matter components opens up the possibility that 
they may be able to explain the current cosmological 
acceleration and possibly dark matter. That these cor- 
rection components become significant only at specific 
scales is encouraging, as this may provide an explanation 
why classical general relativity is valid on sub-galactic 
scales, but requires the input of addition matter on super- 
galactic scales. However it remains to be seen if it is 
possible for such corrections to meet observational con- 
straints. 

This work can be used in two complementary ways: 
top-down and bottom-up. The former is to apply this 
method to any fundamental theory that can produce a 
WDW like equation on large scales and so characterise 
the new, phenomenological effects of the theory. The 
latter would estimate the types of QG corrections that 
would be necessary to produce certain desirable effects 
(e.g., dark energy) at particular scales. The case of 
loop quantum cosmology illustrates how both approaches 
can be used to mutual benefit. The underlying theory 
was used to calculate the form of the QG corrections 
to the WDW equation. The phenomenological conse- 
quence (eventual domination of the correction terms and 
hence a break down in large scale classical physics) high- 
lighted the importance of modelling lattice refinement 
within cosmological models. 

Considering the phenomenological consequences of a 
theory is the first step to testing it against experimental 
and observational data. We have shown here that for QG 
this principle can be invaluable in guiding our search for 
the full theory. 
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